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ABSTRACT 

We present results of Hubble Space Telescope NICMOS H-band imaging 
of 73 of most luminous (i.e., \og\LiR/ Lq] > 11.4) Infrared Galaxies (LIRGs) 
in the Great Observatories All-sky LIRG Survey (GOALS). This dataset com- 
bines multi-wavelength imaging and spectroscopic data from space (Spitzer, HST, 
GALEX, and Chandra) and ground-based telescopes. In this paper we use the 
high-resolution near-infrared data to recover nuclear structure that is obscured by 
dust at optical wavelengths and measure the evolution in this structure along the 
merger sequence. A large fraction of all galaxies in our sample possess double nu- 
clei (~63%) or show evidence for triple nuclei (~6%). Half of these double nuclei 
are not visible in the HST B-band images due to dust obscuration. The majority 
of interacting LIRGs have remaining merger timescales of 0.3 to 1.3 Gyrs, based 
on the projected nuclear separations and the mass ratio of nuclei. We find that 
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the bulge luminosity surface density LBuige/Rluige increases significantly along 
the merger sequence (primarily due to a decrease of the bulge radius), while the 
bulge luminosity shows a small increase towards late merger stages. No signifi- 
cant increase of the bulge Sersic index is found. LIRGs that show no interaction 
features have on average a significantly larger bulge luminosity, suggesting that 
non merging LIRGs have larger bulge masses than merging LIRGs. This may 
be related to the flux limited nature of the sample and the fact that mergers 
can significantly boost the IR luminosity of otherwise low luminosity galaxies. 
We find that the projected nuclear separation is significantly smaller for ULIRGs 
(median value of 1.2 kpc) than for LIRGs (mean value of 6.7 kpc), suggesting 
that the LIRG phase appears earlier in mergers than the ULIRG phase. 

Subject headings: galaxies: active — galaxies: bulges — galaxies: evolution — 
galaxies: interactions — galaxies: starburst — infrared: galaxies 



Introduction 



Luminous IR Galaxies (LIRGs: i.e., Ljr > lO^^Lf;^) were d iscovered in the 1980s with 
the Infra Red Astronomical Satellite (IRAS; ISoifer et al.l Il984[ ). due to their strong emis- 
sion at 12-100//m, and today are an important class of objects for understanding massive 
galaxy formation and evolution. Follow-up observations with ISO and Spitzer have shown 
that LIRGs comprise a significant fraction (50-70 %) of the cosmic in f rared background and 
dominate the star - formation activity sX z > 1 ([Elbaz et al.l l2002l : iLe Floch et al.l 12005 



Caputi et al. 


2007; 


Bridee et al. 


2007: 



20091). Ultra Luminous IR Galaxies 



(ULIRGs: i .e., Lj^ > lO^^Lp) are typically a mixture of disk galaxy pa i rs, interactions 



1996 



Murphv et al 



dust (ISanders et al 



1991 



or ni e rgers (jJoseph &: Wrightl Il985l : lArmus et al.l Il987l : ISanders et al.l Il988l : IClements et al. 



19961). (U)LIRGs are als o observed to be rich in molecular gas and 



Solomon et al.l 1 19971 ): interaction drives the gas inward, resulting 



in enhanced star formation rates. In the local universe, ULIRGs are relatively rare, but are 
about a thousand tir nes more frequen t at z > 2 as shown by observations of the population 
of sub- mm galaxies ( iBlain et al.l l2002l : I Chapman et al.l |2005| ) . 



Recent discoveries have revealed that the central black hole mass is constrained by 



and closely related t o properties of the ho st galaxy's bulge (e.g. iMarconi fc HuntI 12003 



Haring fc Rixl 12004 : iGebhardt et al.l l2000l ). which suggests a co-evolution of the central 



black hole and the galaxy itself. Given both t^ 
evolye into massive ellipt ical and SO galaxies (e.g. 



e evidence that merging 


LIRGs likely 


Barnes & Hernquist 


1992 




Genze 


et al. 



200ll : iTacconi et al.l 120021 ) and that LIRGs are stellar nurseries (see e.g. ISanders et al 



1991 
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Solomon et al.l 119971 : iMagnelli et al.ll2009t iHowell et al.ll2010[ ). the study of nuclear regions 
of a complete sample of LIRGs spanning all merger and interaction stages allows us to 
investigate the coeval nature of both processes. 

The picture that mergers transform disc galaxies into elliptical gal axies is supported 
by numerical simulations (e.g. iBarned Il990t iKormendy fc Sanders! 119921 ). potentially link- 
ing LIRG activity to a key step in the formation of the galaxy morphologies seen locally. 
Other recent simulations have shown evide nce that this may be more complicated (e.g. 
Bournaud et al.ll2005l : [Robertson et al.ll2006l ). suggesting the importance of different merger- 
mass ratios, gas-fraction and feedback processes due to supernovae or AGN activity. More- 
over, the exact number of merger progenitors is contro yersial; while the vast majority of 



inter acting (U)LIRGs are likely two galaxy systems (e.g 
2OO2I ). there may be also a small fraction of LIRG s that involve t 



Murphy et al.lll996l : IVeilleux et al. 



stantial star- formation between the merging nuclei ( Vaisanen et al.l 



iree galaxies and sub- 



20081). While previous 



James et al. 19991 : iRothberg &: Joseph 



groun d-based observations (e.g. IWright et al.lll990 : 
20041 ) have revealed that a significant fraction of merger remnants seems to exhibit radial 
profiles similar to elliptical galaxies based on their large-scale appearance (at radii of several 
kpc), we are able to probe the first time the stellar light distribution in the central kpc for 
a large sample of nearby mergers. 

Merger models that incorporate hydrodynamics and star formation suggest that dis- 
sipation in mergers produces central st arbursts which l e ad to a steepening of the central 
profile well ab ove the extrapolation of a Ide Vaucouleurs I ( 119481 ) profile fit to the main body 



of the galaxy ( Hernquist 1993; Mihos fc HernquistI 



I994J ). Indeed, recent observations have 



revealed such centra! 



( jRothberg fc Joseph 



"ex t ra light" or "cuspiness" components for some elliptical gala xies 



200' 



particular, iHopkins et al. 



Ferrarese et al.ll2006l : ICote et al.ll2007l : iKormendv et al.ll2009f l. In 
(I2OO8I ) have shown that the observed extra light can be identified 



with the central density excess produced in simulations of gas-rich mergers and that suc h 
components are ubiqui tous in the l ocal c usp elliptical populat i on ([Hopkins et al.l l2009al ). 
based on samples from iLauer et al.l (120071 ) and iKormendy et al. J2009h . However, although 
the idea that almost all the gas falls to the center to fuel a central starburst or AGN might 
be true for ULIRGs, this process may be less prevalent in LIRGs. In principle, a LIRG 
sample spanning all merger stages can be used to test the model of merger- induced cusp- 
building alon g the merger stage sequence. The Great Observatories All-sky LIRG Survey 
(GOALS, see lArmus et al.l 120091 ) provides an excellent sample to compare to these models 
and observations. 



To reveal the detailed structure of the nuclear regions, where dust obscuration may mask 
star clusters, AGN and additional nuclei from optical view, we study 73 of most luminous 
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(U)LIRGs in the GOALS sample using the Hubble Space Telescope (HST) Near Infrared 
Camera and Multi-Object Spectrometer (NICMOS). This high- resolution Near Infra- Red 
(NIR) imaging sample is unique not only in its completeness and sample size, but also in 
the proximity and brightness of the galaxies. Moreover, a wealth of additional information 
about the properties of these galaxies (e.g. mid-IR distribution, mid-IR line diagnostics, 
merger stage) is available, since the GOALS project combines multi-wavelength imaging 
and spectroscopic data from space (Spitzer, HST, GALEX, and Chandra) and ground-based 
telescopes. In particul ar, a detailed multi-wavelength study of the LIRG Markarian 266 



( iMazzarella et al.ll2010[ ) reveals nuclear and galac tic-scale outflows, shock- heated gas and the 



presence of a dual AGN. lDiaz-Santos et al.l ( 120101 ) found that at least half of the MIR emission 



is extended for more than 30% of local LIRGs, as well as an increasing compactness towards 
the final stage of major merger interaction. Furthermore, the relationship be tween the IR 



and UV properties (e.g. the IRX-/3 relation as a function of LIRG properties iHowell et al. 



20101 ) has been investigated. The optical GOALS HST images (ACS/WFC F435W and 
F814W filters) have clearly shown the inner spiral structure, dust lanes, exten ded filamentary 
emission, and star clu sters in the nuclear regions and tidal tails of LIRGs (jVavilkin et al. 



2010l : lKim et al.ll2010f l. 



To understand the manner in which clusters, AGN, and nuclei evolve as a function of 
luminosity and merger stage, it is essential to first unveil the hidden nuclei and to model 
the nuclear structure. This requires both high spatial resolution (~100 pc) and the ability 
to penetrate the dust. Unfortunately, it is almost impossible to reveal the circumnuclear 
regions and to identify the nuclei in LIRGs with observations at optical wavelengths, simply 
because they are extremely dusty. On the other hand, NIR imaging is 1) less affected by 
dust extinction and hence provides an almost unobscured view on the nuclear region of 
galaxies, and 2) best suited to trace the old stellar population since NIR light is less biased 
by young blue stars that can dominate the optical light but contribute only a minor fraction 
by mass to the total stellar population. In particular, detailed information about the nuclear 
structure is crucial to reveal double or multiple nuclei to estimate merging time scales which 
can provide important constraints for the lifetime of the IR luminous phase and associated 
galaxy evolution scenarios. Therefore, HST NICMOS images are perfectly suited to observe 
more directly the true nuclear morphologies by combining the high resolution of the HST 
(0.15" FWHM) which corresponds to 30 - 300 pc over the distance range 40 - 400 Mpc (with 
a median resolution of 106 pc at a median distance of 142 Mpc) at wavelengths where dust 
extinction is reduced by an order of magnitude compared to visual wavelengths. 

Our NICMOS imaging has revealed the presence of double and triple nuclei LIRG sys- 
tems. Fitting the NIR light distribution allows us to extract the structural stellar components 
of the galaxies and to study their evolution along the merger stage sequence. Our NICMOS 
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sample and the observations are d escribed in § |2l The modehng of the stellar components 
using GALFIT (jPeng et al.l |2002[ ) and their parameters are characterized in § I3.2.H and 
the results are presented in § [3] together with the numbers of double and multiple nuclei 
found in our sample. In § H] we discuss the number of nuclei that are obscured by dust 
and its implication for high-redshift studies, merger time scales, the evolution of the cen- 
tral stellar morphology along the merger stage sequence, and the role of AGN activity and 
young stellar populations on nuclear properties. For convenience, we refer to the full sample 
{\og[Liji/ Lq] > 11.4) as LIRGs despite the inclusion of a small number (17 systems, ~23%) 
of ULIRGs {\og[LjR/LQ] > 12.0). 



2. Observations and Data Reduction 

The GOALS NICMOS imaging pro gram targets the nu clear regions of all systems in 



the IRAS Bright Galaxy Sample fRGBS: ISanders et al.ll2003l ) with \og[LiR/LQ] > 11.4. The 



REGS is a complete sample of all extragalactic objects with f^{60fim) > 5.24 Jy, covering the 
entire sky surveyed by IRAS at Galactic latitudes \b\ > 5°. The RBGS contains 203 LIRGs, 
of which 88 have luminosities of log[Ljfi/ Lq] > 11.4, the luminosity at which the local space 
density of LIRGs exceeds that of optically-selected galaxies. These galaxies are the most 
luminous members of the GOALS sample and they are predominantly mergers and strongly 
interacting galaxies. Most of the single spiral galaxies or widely separated, weakly/non- 
interacting pairs are filtered out by this luminosity threshold. The fraction of single spiral 
galaxies and non-interacting pairs in our NICMOS sample is ~30%. The redshift range of 
our sample is 0.01 < z < 0.05 and hence the galaxies are bright and well-resolved due to 
their large angular size. Here we present the first results of 73 (56 LIRGs, 17 ULIRGs) out 
of the 88 most luminous (U)LIRGs with a range in luminosit y of (11.4<log[LfB/ Lf7)] < 12.5) 



representing the high-luminosity end of the GOALS sample ( lArmus et al.ll2009l ) 



New HST images with NICM0S/NIC2 have been obtained using the F160W filter for 44 
LIRGs within a total time of 59 orbits (program 11235, Surace P.I.). These data are combined 
with twenty-nine LIRGs that already had high-quality archival NICMOS data (see Tab. [1] 
for an overview of our sample). The data were collected using camera two (NIC2) with a 
field of view of 19.3" x 19.5", a pixelsize of 0.075", and are dithered to yield a total field area 
of typically 30" x 25". Fifteen of the galaxies in our sample have two apparent optical nuclei 
separated by 15"-40", or inclined disks with extensions greater than 20", and thus required 
two pointings per orbit. Note that NICMOS failed during execution of this program, so that 
not all targets in the complete sample of 88 LIRGs were observed. The remainder are being 
observed with WF3, but these data are not yet collected. The data reduction and calibration 
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has been done using the standard HST pipehne. Additional corrections were made to the 
individual frames to correct for bias offsets and we recombined the images using the STIS 
software. A World Coordinate System (WCS) correction could not be performed directly 
using cataloged star coordinates given the limited Field of View (FOV) of the NICMOS 
camera NIC2 (19.2" x 19.2"). Instead the HST ACS images, which have a far larger FOV 
(202" X 202"), have been calibrated with star position references (2MASS) in a first step 
and we subsequently carefully applied a WCS transformation of our NICMOS images given 
the corrected HST ACS images using corresponding reference points in both images (such 
as stars not catalo ged, bright star cluster knots, or similar features). Furthermore, GALFIT 
( iPeng et al.ll2002l ) is used to fit the F160W light from the underlying old stellar population 
and the AGN (see §. EXUfor details). 



3. Results 

Here we present results of the nuclear stellar properties of 73 LIRG systems. Most of 
these systems are major merger, however, some LIRGs in our sample are isolated undisturbed 
galaxies, suggesting that other mechanisms rather than merger interaction may trigger the 
IR luminosity in those isolated undisturbed LIRGs (e.g. high gas fraction and stellar mass, 
cold gas accretion, or other internal dynamical processes). Therefore, these two sub-classes of 
LIRGs, interacting and isolated galaxies, are treated separately in our studies. Furthermore, 
to study the stellar nuclear properties as a function of merger stage (see § 13. 2p . we have 
applied a classification into different merger stages of our interacting LIRGs as described in 
the following. 

The non-merger /merger separation as well as the merger classification scheme is based 
on the wide-field I- and B-band images. While the HST NICMOS H-band images reveal 
the unobscured nuclei and stellar components, the HST ACS I- and B-band images provide 
a larger FOV for the detection of companion galaxies and are very sensitive for tracing 
tidal tails and other interaction features which are essential to identify the merger stage. 
This provides a distinct advantage over using a simple projected nuclear separation to infer 
the merger stage, which suffers from degeneracies when exploring a large range in stages. 
To establish a merger classification we carefully separated at first our sample into galaxies 
that exhibit interaction features (e.g. tidal tails, bridges) or have nearby companions, and 
those without any interaction features or nearby companions (undisturbed single galaxies), 
excluding the latter from our merger sequence analysis. The mergers are c l assified into si x 



different stages (based on the scheme first proposed by lSurace et al.lll998al : ISuracelll998bl ): 



an example for each merger stage based on the I-band images is shown in Fig. [T] 1 - 
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separate galaxies, but disks symmetric (intact) and no tidal tails, 2 - progenitor galaxies 
distinguishable with disks asymmetric or amorphous and/or tidal tails, 3 - two nuclei in 
common envelope, 4 - double nuclei plus tidal tail, 5 - single or obscured nucleus with long 
prominent tails, 6 - single or obscured nucleus with disturbed central morphology and short 
faint tails. In sum, these stages may be broadly characterized into pre-merger (1), ongoing 
merger (2,3,4), and post-merger LIRGs (5,6). While for merging stages 2-6 at least two 
galaxies are obviously interacting or have recently undergone an interaction ( post-merger), 
this must be not necessarily true for stage 1 as there is no evidence that an interaction 
(major merger) is responsible for the LIRG activity (see discussion below). However, the 
maximum line-of-sight velocity difference of our class 1 mergers is very small (ranging from 
65-160 km s~^), which implies that they are very likely bound together as merger. 



3.1. Double and Multiple Nuclei 

We find double as well as multiple nuclei with a large range of nuclear separations (100 pc 
to >10 kpc). Fig. |2] shows three typical merger examples of our NICMOS images including 
a comparison with the mid-IR emission as observed with 24/im MIPS images (Figures for all 
LIRGs are available in the online version of the Journal). We find that a significant fraction 
of our sample shows sign for late merger stage which is typically characterized by tidal tails 
(see §. 13.2.31 for more details), and in some cases, by a disturbed center (with eventually 
unresolved double nuclei). 
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Fig. 1. — The figure at the top shows examples of I-band images representing the classifica- 
tion in different merger stages, ranging from 1 (top left) to 6 (bottom right). The stages (see 
text for details) can be broadly characterized into prc-merging (class 1), merging (2,3, and 
4), and post-merging galaxies (5 and 6). Note that most of our class 1 objects have smaller 
projected separations than this example (AM07020-601), however, a line-of velocity differ- 
ence of 63 km implies that these two galaxies are very likely bound together as merger. 
The figure at the bottom shows two examples of undisturbed isolated galaxies. Mergers and 
isolated galaxies are treated separately thoughout all our studies. 
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Fig. 2. — Three typical merger examples (top: ESO255-IG007, center: ESO060-IG016, bot- 
tom: IRAS 16399-0937) which also show a significant difference in the distribution between 
the NIR light (left: HST NICMOS F160W images scaled with the square root to reveal also 
the outer extent of the galaxies and with contours to identify double and multiple nuclei) 
and the mid-IR hght (right: HST NICMOS image scaled linear and overlaid with MIPS 
24 //m contours) . Figures for the rest of the sample are available in the online version of the 
Journal. 
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To understand how star clusters, AGN, and nuclei evolve as a function of luminosity 
and merger stage, it is essential to first accurately count and measure the central structures. 
Although NIR images trace the old stellar nuclear structure, the possibility of NIR emission 
from accompanying intense star-forming regions rather from a secondary nuclei cannot be 
ruled out (see §. 14.51 for more detailed discussion). We carefully checked different contour 
levels for each galaxy to identify the nuclear structure and possible associated galaxy disks. 
Figures with HST NICMOS contour levels for all LIRGs are available in the online version 
of the Journal. To estimate the number of systems with at least two nuclei we followed a 
conservative approach. We count all LIRGs in our sample that have primarily two galaxies 
with each of them exhibiting a (eventually disturbed or merging) disk as well as LIRGs where 
no clearly visible disk could be associated with every nucleus, but exhibiting two significant 
separated NIR light concentrations (with spherical/elliptical shape expected from nuclei) 
with NIR luminosity ratios < 3 : 1. 

To estimate the number of galaxy pairs that are not within the NICMOS Field of 
View (FOV), we take into account galaxy pairs that are visible in the optical HST ACS 
images (F814W filter, I band) that show signs for interaction such as tidal tails. Out of the 
73 systems, five were interacting systems with one of the companions outside the NICMOS 
FOV. We also find 5 LIRGs with galaxy pairs within 100 kpc that show no strong interaction 
features, and hence are not considered as interacting LIRG systems in our analysis. Note that 
up to 10% of galaxies in our sample that are not classified as double nuclei systems show 
a (slightly) disturbed center (likely caused by a minor merger or strong secular evolution 
effects such as bars or nuclear spirals) but not two separable NIR light concentrations that 
are typical for double nuclei. 

We find that the fraction of LIRG systems (logfL/^j/LQ] = 11.4 — 12.5 ) with at least 
two interacting nuclei is 63% (see Fig. E]), with a median projected separation (see Fig. H]) of 
7.2 kpc. If we also take into account apparent single nuclei in galaxies that exhibit tidal tails 
(given the spatial resolution of 0.15" as upper limit of the projected nuclear separation), the 
median projected separation is 4.9 kpc. Six double nuclei that are visible in the H-band can 
not be identified as two separated nuclei in comparable resolution HST I-band images (ACS 
F814W filter), likely due to dust obscuration (not visible B-band images as well). Triple 
nuclei are found for ESO 255-IG007, II ZW 096, NGC 6670, and possibly UGC 02369 and 
NGC 3690. We note that the fraction of triple nuclei systems in our sample (~6%) is similar 
to the percentage of triple syste ms (~5%) found in a population of ~11,000 galaxies in the 



Local Supercluster and vicinity (IMakarov fc KarachentsevI 120091 ). 



To obtain the number of nuclei that are obscured by dust, we compare the NICMOS 
images with comparable resolution HST B-band images (ACS F435W filter) and find that 
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only ~50% of the double nuclei are visible in the B-band (see Fig. |3]). In detail, the fraction 
of LIRG systems with at least two nuclei visible in the B-band within the NICMOS FOV 
is 28%. Excluding the NICMOS FOV limit we obtain a fraction of 35% LIRG systems 
with interacting double nuclei visible in the B-band. Interestingly, almost half of the dust- 
obscured nuclei (10 out of the 20 that are seen in H-band but not in B-band) are responsible 
for the bulk of the mid-IR emission (MIPS-24/im). A typical example is ESO 099-G004 which 
shows two nuclei in the H-band with a projected nuclear separation of ~3.6 kpc, neither of 
which can be seen in the B-band (see Fig. Interestingly, all of the mid-IR emission is 
emitted from one nucleus in this double nuclei system. Another example is IRAS 16339- 
0937, which shows clearly two nuclei (separation distance ~3 kpc) of roughly the same size 
and luminosity in the H-band, albeit most of the mid-IR emission (>90%) originates only 
from one of the nuclei. The true nature of this source is not evident in the B band data (see 

Fig. ED- 
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□ Apparent single nuclei 



□ Two nuclei but no 
interaction features 

□ Two nuclei (but with one 
outside the NICMOS FOV) 

■ Two or more nuclei 



H-Band l-Band B-Band 



Fig. 3. — Comparison of the number of LIRG with multiple nuclei detected in the H, I and B 
bands. From bottom to top: (1) The number of systems with at least two nuclei in galaxies 
that show interaction features (in dark grey), (2) with two galaxies (visible in I and B band) 
but only one of them within the NICMOS FOV (in grey; we assume that galaxies visible in 
the I band would also been seen in the H band), (3) with double nuclei whose host galaxies 
show no sign for interaction (neither in the H, I, or B band images; in light grey), and (4) 
galaxies with apparent single nuclei {AR < 150pc, in white) and major merger interaction 
features . These results show that half of the double nuclei are obscured by dust as they are 
not visible in the B band. Triple nuclei are found in at least 3 (possibly 5) LIRG systems 
(-6%). 
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Fig. 4. — Histogram of the projected nuclear separation between double nuclei. Fourteen 
apparent single nuclei that exhibit prominent interaction features, such as tidal tails, are also 
included in white, converting our angular resolution limit of 0.15" to the linear size at their 
corresponding distance. For triple nuclei systems, the smallest projected nuclear separation 
is chosen. 



15 25 00'' 24"' 59"' 58=' 57'' 
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Fig. 5. — A three-color composite image (red: H-band, green: I-band, blue: B-band) with 
mid-IR (24 fim) light emission (contours) for ESO 099-G004. This example clearly shows 
the importance of NIR images to reveal multiple nuclei that are not visible in B-band images 
(blue) due to dust extinction although dominating the mid-IR light emission (contours). 
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Right Ascension 
Center; R.A, 16 42 40.21 Dec -09 43 13.5 



Fig. 6. — IRAS 16399-0937; three-color composite image (red: H-band, green: I-band, blue: 
B-band) with MIPS 24 /im contours. This example shows two nuclei (separation distance 
~3 kpc) of roughly the same size and luminosity in the H-band, albeit most of the mid-IR 
emission (>90%) originates only from one of the nuclei. 
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3.2. GALFIT Results 

3.2.1. GALFIT Analysis 
To derive detailed information about the structural properties of LIRGS we pe rformed a 



2-dimensional decompositi on of all galaxies in our NICMOS sample using GALFIT (jPeng et al 



2OO2I : IWeinzirl et al.ll2009[ ). As most of the galaxies in our sample are mergers, and hence 
exhibit disturbed outer disks, we have not attempted to obtain a precise model for the entire 
merger system as arbitrary fitting of the high spatial frequency components does not add to 
a parametric model of the systems. Instead, we focused on an accurate measurement of the 
bulge luminosity and central concentration indices, as these parameters provide important 
information about the build-up of stellar mass in the center of galaxies during the merger 
process. We performed a careful analysis of each galaxy, also identifying cases where a global 
bulge-disk decomposition would not work. 

Running GALFIT on an image requires initial preparation. The desired fitting region 
and sky background must be provided, and the PSF image, bad pixel mask (if needed), 
and pixel noise map must be generated. Our analysis steps are described as follows: (1) we 
converted the integrated galaxy counts per second to an apparent Vega magnitude using the 
NICMOS data calibration description (Dickinson et al. 2002): 

m = -2.5 log{FcPs x PHOTFNU/ F^{Vega)), (1) 

with the source flux Fcps in counts per second, the Vega flux F^[VEGA]) (for NIC2 at 
F160W, Fy[Vega] is 1043.6 Jy), and PHOTFNU as conversion factor from counts per second 
to Jansky (1.49585 x 10^^ Jy s count"^). (2) The fitting region must be large enough to 
include the outer galaxy disk or entire merger system, as well as some sky region. (3) 
Sky backgrounds were measured separately using a region where no emission is present 
and designated as fixed parameters for all GALFIT procedures. In general a fitting of the 
sky background is not recommended as the wings of the bulge Sersic profile (jSersic 1968 ) 



can become inappropriately extended, resulting in a Sersi c index that is too high. (4) We 



simulated the required PSF image with the TinyTim code (iKrist &: Hooklll997l ) using a pixel 
scale identical to that of the observations. (5) The uncertainty images are generated by the 
HST reduction pipeline and used as weighting images required for GALFIT to perform the 
minimization. 

After these preparation steps we carried out an iterative process to perform a 1-3 com- 
ponent decomposition of the NICMOS images described in the following: GALFIT requires 
initial guesses for each component it fits and uses a Levenberg-Marquardt downhillgradient 
algorithm to determine the minimum based on the input guesses. At first we fitted single 
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Sersic profiles that are centered on the galaxies nuclei (typically 1-2 nuclei per image). For 
merger systems that consists of two or more galaxies, we take advantage of GALFITs ability 
to fit simultaneously multiple components at different spatial regions. The Sersic profile has 
the following form: 

S = S,e-^[(^/«^)''"^-^l (2) 

with the effective radius i?e, the surface brightness Sg at Re, the power law index n^, and 
k coupled to Ug such that half of the total flux is always within Rg. The advantage of the 
Sersic profile is that it is very useful for modeling different components of a galaxy such 
as bulges, flat disks, and bars, simply by varying the exponent. As the index becomes 
smaller, the the core flattens faster for R < Rg, and the intensity drops more steeply for 
R > Re- Classical bulges typically have a Sersic index of = 4 (in general known as de 
Vaucouleurs proflle), but can vary between 1 < < 10, while the outer disk follows an 
exponential proflle (n = 1). 

The GALFIT output parameters based on a single Sersic proflle are used as initial 
parameters for the multi-component fltting and are applied as a residual image in the cases 
that more than one galaxy is present per image. The center coordinates have been flxed 
for the subsequent fltting steps. A multi component flt is performed on each individual 
galaxy. In general a double Sersic component is used, while possible additional galaxies 
are subtracted using their single Sersic parameters from step 1. In cases where the galaxy 
shows no signiflcant disturbances we fltted the sum of an exponential and Sersic component 
with flxed position angle. An additional Se rsic component has been fltted when a bar was 



present. As suggested by iPeng et al.l (120021 ). the bar can be well described by an elongated, 
low index Sersic (n < 1) proflle. For each galaxy a variety of different multi- component 
flts are performed (typically 4-7 different flt setups per galaxy) and compared afterwards. 
Fitted models were selected only as long as the model parameters were all well behaved, 
i.e. requiring that all parameters such as bulge radius and Sersic index converged within 
the range of reasonable values (e.g. < 10). The best flt was chosen based on several 
criteria: First, GALFIT calculates the for each model which generally declines as more 
free parameters or Sersic components are given. However, a small does not necessarily 
mean that the solution with an extra bar component is more correct physically. Thus, an 
increasing for a model with an extra bar component was interpreted that the flt should 
not be adopted, but a decreasing was not considered as a sufficient condition to adopt an 
extra bar component. In cases with prominent bars, a symmetric light distribution due to 
unsubtracted bar light was often found in the bulge-disk residuals. For these cases we added 
an extra bar component since it would otherwise lead to an extended bulge. 

For objects with central point sources such as an AGN, the bulge Sersic index can be 
overestimated unless an extra nuclear component (PSF component) is added to the model. 
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On the other hand, nuclear stellar cusps can mimic a central point source and adding a 
nuclear PSF would lead to an underestimation of otherwise a steeper Sersic index. Therefore, 
we used GALFIT in a first attempt without a central PSF component and calculated the 
central residual excess light from the residual image of the final fit. In cases where a central 
residual was present (for ~10% of the sample) we reran GALFIT with an additional PSF 
component to take into account the extra central light originating from a central point 
source. We also checked all LIRGs that show signs for AGN activity (based on 6.2 /im PAH 
Equivalent Width, see § 14. 5p whether a PSF component is present or not. For those cases 
where the model successfully converged with the extra component, the images were visually 
inspected to verify the presence of a central bright source. Where new model parameters 
were not unreasonable and not identical to the case without the point source, the new model 
was adopted. This was the case for 14 galaxies (17% of the sample), but only for 9 of 
them reasonable bulge parameters (with L^uige, -Reuige, and Ug > 3a) could be obtained. 
For consistency we also checked galaxies with less significant residual light and larger PAH 
EQWs (no significant AGN contribution), but an extra nuclear component was not found 
for any of those galaxies. 



3.2.2. Derived Bulge Parameters and Black Hole Mass Estimation 

The following morphological parameters are derived via fitting of our H-band images 
using GALFIT (see § 13.2. ip : (1) the total bulge luminosity i^Buige, (2) the Sersic index Us of 
the bulge component radial profile, and (3) the bulge radius -Rsuige (see Tab. [2]). We have 
obtained bulge parameters for 73 nuclei (with -Leuige, -Rsuige, and rig > 3(t) i n our sample and 



calcu lated the central black hole masses using the Marconi-Hunt relation (IMarconi fc Hunt 



20031 ) given as 

logMBH = 8.19(±0.07) + 1.16(±0.12) ■ [\og{LH,Bulge) - 10.8] (3) 

with the bulge luminosity log{LH,Buige) measured from the H-band. Moreover, we calculated 
the bulge luminosity surface density defined as LBuigeZ-Rluige' ratio of the core excess 

luminosity to bulge luminosity i^excessZ-^Buige (measured within the central ~ 30 — 300 pc ). 
The central residual light or core excess luminosity /^excess is derived from the center of the 
residual image (in detail the bulge model subtracted H-band image). In case of a significant 
central residual light we subsequently refined the core excess luminosities by fitting a central 
PSF component in GALFIT (see §. 13.2.11) . Two simple examples of GALFIT models and 
residual images are shown in Fig. [7] (GALFIT model and residual images for all LIRGs are 
available in the online version of the Journal). 

We find only a small trend (~ 3o" significance level, see Fig. [8]) in which the bulge 
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Fig. 7.— GALFIT output images. Top left: Region of NICMOS image that is fitted (box 
correspond to fitted region), Top right: Galfit model, Bottom left: Residual image (shows the 
difference between model and NICMOS image) with the same brightness level as the NIC- 
MOS image, Bottom right: Residual image scaled to its maximum brightness (the brightness 
ratio of NICMOS image to scaled residual map is shown as well at the bottom) to highlight 
the faint diffuse emission, spiral structure, and clusters remaining in the model-subtracted 
data. The axes are in scales of arcsec and the figures are orientated in the observational 
frame. Figures for the rest of the sample are available in the online version of the Journal. 
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lumino sity is correlated wit h the total IR luminosity (using the flux densities reported in the 
RBGS. ISanders et al.ll2003l ) as given by the following flt parameters: 



\ogiLBuige/Le) = 3.84(±2.14) + 0.61(±0.18) ■ log(L,^/Le 



0; 



(4) 



This result indicates a small dependence between bulge stellar mass and total IR luminosity 
since a possible contribution of AGN activity on the bulge NIR luminosity is unlikely (as 
discussed in § 14.51) . 
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Fig. 8. — Bulge luminosity deri ved from GALFIT v ersus total IR luminosity (using the flux 
densities reported in the RBGS. ISanders et al.ll2003[ l. Bulges that belong to the same merger 
system are highlighted with an outer ring at a common IR luminosity. A linear flt over all 
data points (dashed line) reveals log{L Bulge /Lq) = 3.84(±2.14) + 0.61(±0.18) ■ logiLm/LQ). 



3.2.3. Bulge Properties as a Function of Merger Stage 

To investigate the effect of the merger process on the central evolution of galaxies, 
such as the growth of the central black hole (BH) or stellar concentration, we study several 
morphological parameters as a function of merger stage: BH masses and bulge luminosities 
(Fig. [9]), bulge Sersic index (Fig. [10]), bulge radius (Fig. [11]), and the bulge luminosity surface 
density (Fig. [T2]) . An overview of the mean, median, mean systematical error (from GAL- 
FIT), and the standard deviation of the bulge properties for all merger stages is presented 
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in Table [31 To compare the bulge properties between different merger stage populations 
we performed Kolmogorov-Smirnov (KS) and Mann- Whitney U (MWU) tests (presented 
in Table H]) which are based on the Cumulative Distribution Functions (CDF) of the bulge 
properties for each merger stage (see Fig. [T31) . 

As shown in Fig. |9l the BH masses (bulge luminosities) show an increase from merger 
stage 3 to merger stage 4-6 by a factor of 1.5-1.8 in their mean and median values. The 
KS and MWU test results confirm that the distributions of the values in each of these bins 
differ from each other (at a probability level of ~90%, see results in Table S]). The bulge 
Sersic index (see Fig. [TOj) tends to decrease from merger stage 1 (mean Ug = 4.8) to 5 (mean 
Us = 2.2). However, a significant change of the distribution of the bulge Sersic index along 
the merger stage populations 2-6 is not supported by the KS and MWU tests. We find 
that the bulge radius Rsuige decreases significantly along all merger stages (shown in Fig. [TT] 
and Table H]) from a median radius of 2.2 (merger stage 1) over 0.8 kpc (merger stage 2-3) 
to 0.3 kpc (merger stage 4-6). With decreasing bulge radius, the bulge luminosity surface 
density L Bulge /R'suige increases about a factor of 5-60 (median values) along the merger 
sequence (1-6, except 3) as shown in Fig [Inland statistically verified in Table ID 

Interestingly, the LIRGs in our sample that show no interaction features (single isolated 
galaxies and merger stage 1 galaxies), have on average a significantly larger bulge luminosity 
(by a factor of ~2), bulge radius (a factor of ~1.5), bulge Sersic index (a factor of ~1.5), 
and a smaller luminosity surface density (a factor of ~10) than LIRGs in merger stage 2-6. 
This is also confirmed by the KS and MWU test results and we will discuss possible reasons 
in §1131 

We find that ~20% of our sample exhibit a core excess to total bulge luminosity 
-^excess /-^^Buige > 0.01, but there is no clear correlation between /^excess /-^Buige and merger 
stage (see Fig. [H]). We find that two non-interacting LIRGs (AM 0702-601 North and 
IRAS 19542+1110) have a very large core excess to total bulge luminosity -^v excess /-Z^Buige 
which is seen as a very bright central cusp in the center of these galaxies. A possible contri- 
bution of AGN activity to the core excess light is discussed in § 14.51 
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Fig. 9. — The distribution (cross marker) of the central black hole masses (left label) along 
the merger stage sequence. The black hole masses are deriv e d from the H-band bulge lumi- 
nosities (right label) using the relation of iMarconi &: HuntI ( 120031 ) . The median and mean 
values of the central black hole mass for each bin are marked with a diamond and circle, 
respectively. The stages (see text for details of the merger classification scheme) can be 
broadly characterized into pre-merging (1 in white), merging (2,3, and 4 in light grey), and 
post-merging galaxies (5 and 6 in dark grey; the nuclei have merged but the structure in 
the disk and tidal tails indicate the source has gone through a merger). Note that we have 
combined merger class 5 and 6 as post-merger stage LIRGs. 
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Fig. 10. — The bulge Sersic index Ug along the merger stage sequence. The Sersic index 
defines the degree of curvature of the radial brightness profile (see Equation [2]). A smaller 
value of Us indicates a less centrally concentrated profile and a shallower (steeper) logarithmic 
slope at small (large) radii. Merger classification and markers are defined as in Fig. [91 
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Fig. 11. — The bulge radius Rsuige along the merger stage sequence. Merger classification 
and markers are defined as in Fig. |9l 



- 25 - 



+ 
+ 

♦ 

+ 

+ 



N ^ 6 



+ 

i 



+ 
+ 
+ 



N ^ ?? 



+ 

t 

+ 
+ 

+ 

N ^117 



+ 
+ 

I 



N ^ 10 



+ 

+ 
+ 

f 
4 



+ 

+ 



N ^ 1.^ 



3 4 
Merger Stage 



Fig. 12. — The bulge luminosity surface density as defined as L Bulge / ^^Buige along the merger 
stage sequence. Merger classification and markers are defined as in Fig. |9l 
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Fig. 13. — The Cumulative Distribution Functions (CDF) of the bulge luminosity (top left), 
effective bulge radius (top right), bulge Sersic index (bottom left), and bulge luminosity 
surface density (bottom right) for each merger stage (legend; see text for detailed description 
of merger stage classification). Note that we have combined merger class 5 and 6 as post- 
merger stage LIRGs. 
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Fig. 14. — The core excess luminosity fraction defined as the ratio of the core excess luminos- 
ity to bulge luminosity L excess / L Bulge along the merger stage sequence. Merger classification 
and markers are defined as in Fig. |9l 
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3.3. Non-Nuclear Mid-IR emission 



Although most of the LIRGs have peaks in their 24/im MIPS emission that correspond 
to peaks in their H-band NICMOS emission (i.e. the stellar nuclei), there are a few LIRGs 
with warm dust emission that is clearly offset from the stellar nucleus. A typical example for 
such a galaxy is II-Zw-096 which clearly shows that 80% of its total infrared luminosity comes 
from an extre mely compact red source not obviously associated with the nuclei of the merging 
galaxies (see llnami et al.l 120101 ). Other similar galaxies in our sample are IC 1623 (>80% 
of MIPS-24yum from a secondary nucleus), ESO 593-IG008 (>80% of 24/im originates from 
possible secondary nucleus), and IRAS F14348-1447 (24/im emission peak between double 
nuclei). One possible explanation is that these galaxies are presumably not yet relaxed with 
off- nuclear starbursts and/or strong shocks. Moreover, three LIRGs show also dominant off- 
nuclear 24/im emission that is very likely associated with extended star formation in spirals 
arms: Arp 256 North (>70% of 24/im in spiral arm), NGC 5257 (>80% from outer spiral 
arm), and NGC 6670 West (~50% from spiral arm). In particular one LIRG double system, 
namely NGC 5257, stands out from the rest of our sample, as it shows that almost all of its 
mid-IR emission is generated at the end of one of the spiral arms (in both galaxies), ~7 kpc 
away from the center (Fig. [151 note that only one galaxy is shown due to the limited FOV 
of NICMOS). This seems to be very surprising given their very regular spiral disks and that 
no visible counterpart is seen at the center of the mid-IR emission, neither in the H- nor 
in the B-band image. This phenomenon might be also associated with Tida l Dwarf Galaxy 
(TDG) formation, w' hich is also suggested for other interacting galaxies (jPuc &: Mirabel 
19981 : IDuc et al.ll2000l ). Since NGC 5257 consists of two galaxies with a projected separation 
distance of ~50 kpc, the exact cause of the enhai iced disk emission is not clear. Similar cases 
are well known for, e . g., the Ant ennae galaxies (IMirabel et al.lll998l ) and the LIRG Arp299 
Jcharmandaris et al.l[2002l . l2004h . 
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Fig. 15. — NGC 5257; three-color composite image (red: H-band, green: I-band, blue: 
B-band) with MIPS 24 fim contours. 
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Discussion 



4.1. Obscured Nuclei 



Understanding the exact number of merger progenitors for LIRGs is important as 
these numbers provi de information abo ut merger timescales and subsequently LIRG ac- 
tivity timescales (see iMurphy et al.lll996l ). The study of the small and large-scale morpholo- 
gies of LIRGs in the nearby universe allows us to understand the physical processes that 
drive galaxy evolution. While ULIRGs are predominantly merger systems, the fraction of 
merger systems and the m erger timescales for L I RGs are more controversial: The imaging 
analysis of 30 local LIRGs lAlonso-Herrero et al.l ( 120061 ) indicates that most of these LIRGs 
(\og[Liji/LQ]= 11.0 — 11.9) have prominent spiral patt erns, and of those, a non negligible 
number are weakly interacting or even isolated systems. iBridge et al.l (120071 ) find in an anal- 
ysis of MIPS 24 yum detected and undetected mergers (0.2 < -2 < 1.3), that a larger fraction 
of LIRGs appear later, in the merger phase, than in the pre-merger, close pair phase. How- 
ever, it is expected that the fraction of major mergers increases with IR luminosity. Thus 
far, the trigger mechanisms and end products of LIRGs have not been constrained as well 
as those of local ULIRGs. 

We find for our sample that the fraction of LIRG systems (log[L/ij/L0]= 11.4 — 12.5) 
with at least two interacting nuclei is 63% (see Fig. |3]). Note that the fraction of major 
merger for the entire GOALS sample (\og[Liji/LQ]> 11) is ~50% (based on Spitzer IRAC 
and HST data). Interestingly, the comparison with our HST ACS B-band images revealed 
that roughly half of the double nuclei seen in the HST NICMOS H-band are obscured by dust 
(not visible in the B-band). Therefore, NIR observations of LIRGs and ULIRGs at high- 
redshift (z>2), which correspond to the B-band rest-frame range, might require significant 
correction factors of ~2 or higher to accurately estimate the true number of multiple nuclei 
systems. 



In particular, IPasvra et all (l2008l . using HST NICMOS F160W) found that the fraction 
of binary systems in ULIRGs at redshift 2 is roughly 50%, which is significantly smaller than 
that of the local ULIRG populati on (nearly 100% of objects in the IRA S 1 Jy sample show 
signs of rnajor rn erger interaction iKim et al.l l2002l : IVeilleux et al.l 120021). Also in studies of 



Bell et al.l (l2005l . MIPS 24 /xm and HST ACS F850LP filter) and Melbourne et all (12005 



MIPS 24 ij,m and HST ACS Viz band) major merger appear to make up a smaller fraction 
of z~ 0.7 LIRGs compared to local systems. The strong wavelength dependence on the 
number of visible nuclei may explain some of the apparent discrepancy between the (U)LIRG 
population at local and high-redshift, particularly in identifying late stage mergers when tidal 
tails have faded and double nuclei are surrounded by dust. 
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4.2. Merger Time Scales 



The space density of LIRGs combined with an estimate of the time-scales of their power 
sources (star formation and/or AGNs) provides important constraints for galaxy evolution 
scenarios that lin k the IR-luminous p hase with the luminous, unobscured AGNs observed 
in Quasar Hosts (ISanders et al.lll988l ). The merger time scale depends on the separation 
be tween the mult i ple nu clei observed in our NIR images (Fig. H]) and we follow the approach 
by iMurphy et al.l (119961 ). namely: In a first stage we adopt a constant velocity Vr until the 
nuclei are 10 kpc apart. Because of the small range in nuclear separations, and the large 
uncertainties involved, we treat the radial velocity as constant. In this case, the time from 
observation until the nuclei are within 10 kpc is (r — 10 kpc)/fr with r as the projected 
nuclear separation. The time spent in the second stage is estimated by calculating the 
dynamical time scale times the mass ratio of the two nuclei, given as (M1/M2) ■ li^r jvc with 
the mass ratio M^/M^, the projected nuclear separation r and the orbital velocity Vc (see 
Murphy et al.l Il996l : iBinney fc Tremaind |2008| ) . As the nu clei spiral together, the re lative 
velocity is similar to the circular velocity at that radius ( Binney fc Tremaine 20081 ). In 
the calculations that follow we adopt an average radial velocity Vr and an a verage circular 



veloc ity fc, of 200 km s ^ for all the double nucleus galaxies in our sample (IMurphv et al. 



1996 ). These average velocities are supported by simulations (e.g. see Fig. 15 in IXue et al. 



20081 ) and the distribution of the line-of-sight velocities differences between the galaxy pairs 
in our sample (typically between 60-180 km s~^; note that the radial velocity is on average 
a factor of 1.5 larger than the line-of-sight velocity component). It is important to note that 
for an individual galaxy all we can estimate from the images and these simple dynamical 
arguments is the time from the present until the time of final merger, when the nuclei of 
interacting galaxies coalesce. We must rely on the distribution of these time scales within the 
sample, along with the morphologies, to allow a calculation of the average remaining merger 
time scales of our sample. For our LIRG sample (log[L/ij/L0]> 11.4) we find a median 
remaining time until the nuclei merge (t — tm) of 4.3 ■ 10* yrs. 

As shown in Fig. [111 we find two peaks in the numbers of LIRGs as function of the 
remaining merger time scale: The largest occurs at 0.3 Gyr< [t — tmerg] <1.3 Gyr, roughly 
representing the first passage of interacting galaxies, while the second peak occurs at the final 
coalescence of the nuclei (most at [t — tmerg] < 10'' yr), including all galaxies with apparent 
single nuclei that show interaction features. Note that the second peak represents the same 
number of LIRGs a nd ULlRGs (each 6). The fact that w e see two peaks is consistent with 
merger simulations (jCox et al.l l2008l : iHopkins et al.l l2008l ) . which show two star formation 
rate (SFR) peaks as function of merger time. Furthermore, the time-scales of these two 
peaks are roughly consistent with the two SFR peaks as seen the in merger simulations, but 
we see slightly more LIRGs at [t — tmerg] < 1 Gyr. In detail, the peak at [t — tmerg] ~ 1 Gyr is 
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Fig. 16. — Histogram of the remaining merger time scale, t — tmerg (see text for details) with 
a median time scale of 4.2 ■ 10® yrs. Fourteen apparent single nuclei in galaxies that exhibit 
prominent interaction features, such as tidal tails, are also included in white, converting 
our angular resolution limit of 0.15" to the linear size at their corresponding distance and 
assuming a nuclei mass ratio of 1:1. 



a bit broader (0.3 Gyr< [t — tmerg] <l-3 Gyr) than we would expect from the merger models 
(0.8 Gyr< [t—tmerg] <l-3 Gyr), indicating that we observe more LIRGs at shorter timescales. 
The fraction of interacting LIRGs is larger for the peak at ~1 Gyr (53% of interacting LIRGs) 
than for the peak at the nuclear coalescence (26%), which likely corresponds to the broader 
width (~0.7 Gyr) of th e SFR peak at ' -- '1 Gyr than the nar row peak (~0.3 Gyr) at the 



nuclear coalescence (see ICox et al.l 120081 : iHopkins et al.ll2008l ). As our GOALS sample is 



a flux-limited snapshot of galaxies passing through the LIRG and ULIRG phase, the data 
suggest that galaxies spend a longer time in the LIRG phase, which is more easily associated 
with the first peak in the model SFR at [t — tmerg] < 1 Gyr. Furthermore, the IR luminosity 
is slightly larger for LIRGs at [t — tmerg] — (median \og[Liii/ Lq] = 11.7) than for LIRGs at 
[t — tmerg] — 1 (median \og[LiR/LQ] = 11.9). This is consistent with our finding that LIRGs 
in late stages of merging have significantly larger total IR luminosities (roughly a factor of 
two) than pre- or non-merging LIRGs. We also discuss this point below (see Fig. [T7]as well). 
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To estimate the timescale that LIRGs spend in the final nuclear coalescence (post-merger 
stage), we have applied the following approach: In principle, the post-merger time scale can 
be estimated by multiplying the merger time-scale with the ratio of galaxies with single nuclei 
to double nuclei. However, this approach might be oversimplified since not all LIRGs might 
have undergone a merging process, and instead other mechanisms may be responsible for a 
significant fraction of LIRGs. In fact, 23% of LIRGs in our HST NICMOS sample show no 
major interaction features in the H-, I-, or B-band (merger stage and 1, see § l3.2p and hence 
it is unlikely that they have undergone a previous major merger. Therefore, we take into 
account only galaxies that are classified as mergers (merger stage 2 — 6, see § 13.21) . The time 



spent in the pre-merger stage for a random galaxy is statisti cally roughly twice t he average 



observed time required to complete the merging t — tm (see iMurphy et al.lll996l ). or about 
lO^yrs for our sample. Thus, the time spent in the post-merger stage is the median merger 
time of LIRGs with double nuclei (2 x 0.43 Gyr) times the ratio of the number of single 
nuclei with interaction features (15 LIRGs) to double nuclei (45 LIRGs). This results in an 
average post-merger time of ~ 2.9 • 10^ yrs, which is very si milar to the typical timescale of 



the S F burst of 0.3 Gyr as seen in major merger simulations (jCox et al.ll2008l : iHopkins et al. 



20081). 



Since we expect most ULIRGs to pass through a LIRG phase, we might expect to see 



a correlation of me rger stage, nuclear separation and IR luminosity (IMurphy et al.l 12001 



Veilleux et al.l 120061 ) . A comparison of LIRGs and ULIRGs in our HST NICMOS sample, 
presented in Fig. [TTl reveals that the projected nuclear separation is significantly differ- 
ent between both populations (probability level of 98% in KS test) with a mean (median) 
projected separation of 3.4 kpc (1.2 kpc) and 11.1 kpc (6.7 kpc) for ULIRGs and LIRGs, 
respectively. 

The results for our subsample of ULIRGs {\or\ LtrI Lf:^] > 12.0) are also in good agree- 
ment with a study for a larger sample of ULIRGs (IMurphy et al.l 119961 ). which is based on 
ground-based observations (seeing of AR ~ 0.8" for K band images), and hence the lower 
spatial resolution limits the ability to resolve double nuclei. Including galaxies with apparent 
single nuc lei (upper limi t s on t he separation of double nuclei), most of the ULIRGs in the 
sample of iMurphy et al. Jl996h ( more than 60%) have nuclear separation between 1-4 kpc 
(median of ~1.8 kpc) which is very similar to our results for the ULIRG popu lation (median 
of ~ 1.2 kpc). Note that we have an overlap of ten ULIRGs with the sample of iMurphy et al. 



(119961 ) and we found for one ULIRG, namely IRAS F19297-0406, two nuclei with a projected 
nu clear separation of 0. 77" that was previously characterized as an apparent single nucleus 
by IMurphy et al.l (119961 ) due to limited spatial resolution of the R-band image (1.9"). How- 
ever, only our entire sample covering a IR luminosity range of log[Lifi/LQ] > 11.4, allows us 
to see a strong decrease of the projected nuclear separation as a function of IR luminosity. 



as shown in Fig. [T71 
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Fig. 17. — The projected nuclear separation as a function of IR luminosity for the observed 
double nuclei (plus marker), as well as apparent single nuclei (arrows down, converting our 
angular resolution limit of 0.15" to the linear size at their corresponding distance) whose 
host galaxies exhibit interaction features. The mean and median values of the projected 
nuclear separation are shown as solid and dashed horizontal lines, respectively. 
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4.3. The Evolution of the Central Stellar Structure along the Merger Stage 

Evidence suggests that ULIRGs may evolve into elli ptical galaxies once the starburst 



subsides and the gas is either used up or expell ed in a wind (ISanders et al.lllQSSt iHeckman et al 



19901 : iGenzel et al.ll2001t IVeilleux et al.ll2006l ). Therefore, the stellar surface brightness pro- 
files of LIRGs should provide a glimpse of the process of bulge and black hole building. Un- 
like previous studies, which focused primarily on ULIRGs {\og[LiR/ Lq] > 12.0), the GOALS 
NICMOS sampl e targets the nuclear regions of all merger systems in the IRAS Bright Galaxy 



Sample (RGBS; ISanders et al.l 120031 ) with \og[LiR/ Lq] > 11.4. Because of the low red-shift 
range of our sample (0.01 < z < 0.05), the galaxies are bright and well-resolved due to their 
large angular size. 

Although we find a significant increase in the bulge luminosity surface density (a factor 
of 5 to 60) along the merger sequence, the growth of the BH mass (bulge luminosity) toward 
later merger stages is only marginal: A factor of ~ 1.8 in bulge luminosity from merger stage 
3 to 5/6 with a probability level of ~90% that both merger stage populations are not drawn 
from the same parent population (KS and MWU test, see Table H]). Two possibilities may 
explain why we do not see a a more significant increase of the bulge luminosity: First, as 
the merging proceeds, the bulge may get partially disrupted or stripped and may assemble 
again after the LIRG phase has passed. Thus, our measured bulge luminosities may miss a 
significant bulge fraction at the late merger stage, resulting in smaller estimated black hole 
masses. A simpler explanation could be that the fiux limited nature of our HST LIRG sample 
(log[L/ii;/L0] > 11.4) excludes galaxies with smaller bulge masses in single isolated galaxies 
and galaxies in the pre-merging stage. Indeed, we find evidence for this, as the LIRGs in our 
sample that show no evidence for interaction features (single isolated galaxies and merger 
class 1) have on average a significantly larger BH mass (a factor of two). Since LIRGs are 
drawn from the top end of the IR luminosity function and mergers boost IR luminosity, we 
might expect the non-interacting LIRGs to be more luminous and hence more massive than 
the progenitors of interacting LIRGs. The interacting LIRGs would tend to include more 
systems with smaller bulges and then grow their bulges during the merging process. If this 
is the case, the non-interacting LIRGs would tend to have more massive bulges - exactly as 
we observed in our sample. A detailed comparison of the bulge masses of non-interacting 
spiral galaxies with log[L/ii;/L0]< 11.4 would be needed to test this hypothesis. 



4.4. Merger-induced Building of Central Starbursts: Comparison to Models 



Dissipation in mergers can generate central starbursts, imprintin g a central "extr a light" 
component into the surface brightness profiles of merger remnants ( iHernquistI Il993l ) . Such 
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an excess of the central light or "cuspiness" in the surface brightness profiles has been found 



for some elliptical ga laxies ( iFerrarese et al.l l2006t ICote et al.l 120071 : iKormendy et al.l 12009 



Hopkins et al.ll2009a( ). This can be explained by the assumption that the envelopes of cusp 
ellipticals are established by violent relaxation in mergers acting on stars present in gas-rich 
progenitor disks, while their cen ters are structured by the relics of dissipational, compact 
starbursts (IHopkins et al.ll2009br). Given the evidence that ULIRG s likely evolve into massive 
elliptical and SO galaxies ( iGenzel et al.ll2001t iTacconi et al.l 120021 ) and the fact that LIRGs 
and ULIRGs host powerful starbursts, the study of our LIRG sample as a function of merger 
stage allows us in principle to test the build-up process of nuclear starbursts. 

To test such a scenario we studied the following parameters as a function of merger 
stage (1 — 6, see §. 13. 2p : (1) the effective bulge radius Rsuige, (2) the bulge luminosity 
L Bulge, (3) the bulge Sersic index Ug which defines the steepness of the inner bulge profile 
(see Fig. [10]), (4) the bulge surface brightness defined as L Bulge / R\uigei ^^cl (5) the central 
light concentration given by the ratio of core excess luminosity to total bulge lu r ninosit y 
(see Fig. [T^ . We have compared our results with recent models of iHopkins et al.l (j2009al ). 
which are based on N-body simulations and smoothed particle hydrodynamics, account for 
radiative cooling and for heating by a UV background, and incorporate a subresolution 
model of a multiphase inters tellar medium ( ISM) t o describe star formation and supernova 
feedback. Since the model of iHopkins et al.l (l2009al ) focuses on local gas-rich mergers (unlike 
high-redshift simulations which have a factor of 2-5 larger gas content), it is ideally suited 
for a comparison to our data-set of local major mergers. 



4.4- i- Bulge Radius and Surface Brightness 



The merger model of iHopkins et al.l (l2009a] ) predicts a decrease of the effective bulge 
radius along the merger process due to gas inflow: Tidal torques excited during major merger 
provide the fuel to power int ense starbursts boostin g the concentration and central phase- 
space density. As shown in iHopkins et al.l (l2009ai ) a gas inflow of e.g ~10% shrinks the 
effective bulge radius to about half its previous size and subsequently enhances the bulge 
surface brightness. Our findings show the same trend as the model predictions (see Fig. [11] 
and Fig. [12]), given a decrease (increase) of the effective bulge radius (surface brightness) by 
a factor of ~ 2 — 3 (5 — 20), and hence provide evidence for the idea that gas inflow can 
dramatically reduce the apparent bulge radius during a merger. 
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4.4- 2. Bulge Sersic Index 



Due to the violent relaxation of stars during a merger, the model of iHopkins et al. 



(120081 ) predicts an increase of the bulge Sersic index Ug towards late merger-stage (from 
t — tmerg = 0.5 Gyr to t — tmerg = Gyr). We do not find a significant increase of Ug (see 
Fig. [TU]and Table H]). This discrepancy can have several reasons, such as that young stellar 
populations and age (metallicity) gradients can shift Ug significantly and/or tidal interactions 
are changing during the more active /merging phases. 



4.5. The Role of AGN Activity and Young Stellar Populations on Nuclear 

Properties 

Throughout most of this paper, we have assumed that the near-infrared H-band data is 
an accurate tracer of the stellar mass, as it is less affected by dust extinction (than the UV or 
optical), and has less of a contribution from young, blue stars. However, LIRGs and ULIRGs 
are known to be powered by starbursts and AGN, which can, in some cases, significantly 
contribute to the near- infrared light. In the following sub-sections we discuss the possible 
contributions of central starbursts and AGN to the near-infrared light. 



4-5.1. Young Stellar Populations 



Evidence suggests that asymptotic giant branch (AGB) and red supergiant (RSG) 



stars can dominate the NIR in the center of so me star burst galaxies (lArmus et al.l Il995 



Rothberg fc Fischerll2010l : [Melbourne et al.ll2010l ). Since high-spatial resolution spectroscopy 
is required to identify and resolve the young stellar populations, their spatial distribution is 
unknown for our sample. However, two scenarios can constrain the possible effects of young 
stellar populations on our measured nuclear properties: (1) If young stars form in the core 
of the galaxies at scales < 300 pc, these populations would be taken into account by our 
measured PSF component and hence not affect our measured bulge properties. (2) In case 
young stars are more widely distributed on scales typical for stellar bulges (median effective 
bulge radius for our sample is 0.7 kpc), they could affect our measured properties. If young 
stars form in a disk, they would manifest an exponential radial profile, and a central GAL- 
FIT component with a Sersic index Ug — 1 (exponential disk) would be expected, instead 
of > 1 (more typical for old stellar bulges), which we find for most of our galaxies (see 
Fig. [TTl) . If this is not the case, namely that the majority of young stars form at bulge 
scales (> 300 pc, not in in the core), are distributed in a spheroid rather than in a disk, 
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and dominate the NIR, then the H-band bulge luminosity may not be a reliable estimate 
of the central BH mass. However, this does not affect our comparison of the central stellar 
structure along the merger stage to merger models (§. l4.4p, since young stellar c ontributions 
are already taken into account in the merger simulations (iHopkins et al.ll2009al Jbl). 



4.5.2. AGN Activity 



As described in § I3.2l we find for some galaxies a large ratio of core excess to total bulge 
luminosity i^excess/-^Buige- One possible explanation of such a bright central light emission 
could be that the LIRG activity is caused by AGN heating. In general, lower Equivalent 
Widths (EQWs) of the 6.2 (j,m PAH feature are associated with AGN activity (see e.g. 
Genzel et al.lll998l : ISturm et al.ll2000t lArmus et al.l 120071 : iDesai et al.ll2007| ) because the hot 
dust continuum in creases and the hard AGN photons may destroy/ionize the PAH molecules. 
Petric et al.l ( I2OIOI ) estimated that AGN are responsible for ~12% of the total bolometric lu- 
minosity of local LIRGs based on several mid-IR line diagnostics measured with the Infrared 
Spectrograph on Spitzer. 

We find no significant trend between the ratio of co re excess to bulge luminosity and 
the PAH EQW (see Fig. [181) of the 6.2 /im PAH feature jPetric et allboiol ). However, one 
example of a galaxy where a central AGN might dominate the NIR emission is AM 0702-601: 
A galaxy in our sample with a very small PAH EQW (0.037 ^m) but also exhibits one of 
the most significant core excess luminosity fractions (i^excess/-^Buigc = 1-5). Counter exam- 
ples, where the core excess light fraction in the H band and PAH EQW are not correlated: 
ESQ 060-IG016 and NGC 3690 East, which have very small PAH EQWs (< 0.27 /im) and 
hence likely large AGN contribution, but they exhibit no significant core excess light fraction 
(LexcessZ-Z^Buige < 0.001). One possiblc explanation might be that those AGN systems are 
so deeply embedded in dust that they don't show a central light excess in the H-band, in 
contrast to som e ULIRGs which show a very large light excess in the center when an AGN is 
present (see e.g. lSurace fc SanderslE999[ ). We find also galaxies (ESO 239-IG002, WKK 2031, 
and NGC 1614) with a significant core light excess (i^exccss/-^^Buigc > 0.05) but relatively large 
PAH EQWs (0.4 - 0.7 /im), suggesting that those galaxies build up a concentrated stellar 
"cusp" in the center due to luminous nuclear starbursts. The sources with the largest core 
excess ( /^excess /-^Buige > 0.05) are not responsible for the decrease in the bulge radius along 
the merger sequence, as the central PSF (even if it is due to an AGN) is already taken into 
account by the GALFIT fitting process (see §. 14. 4p . We obtain roughly the same results for 
the bulge properties along the merger sequence if we exclude the galaxies with small 6.2 /xm 
PAH EQW (EQW < 0.27 /im) from our sample. 
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Lum central residual / total bulge 



Fig. 18. — Equivalent Width (EQW) of the 6.2 /xm PAH feature as function of the core 
excess luminosity to total bulge luminosity. A smaller EQW indicates a larger contribution 
of AGN generated emission to the mid-IR emission. No significant trend is visible between 
6.2 fim PAH EQW and core excess luminosity (see text). 

AGN activity might also effect the observed NIR luminosity and hence one would expect 
a possible contribution to the measured bulge luminosity. This could subsequently explain 
the increase of the BH mass (bulge luminosity) with larger IR luminosity (see Fig. |H] and 
Eq. H]). However, a possible contribution of AGN activity to the bulge luminosity is not very 
likely since most of our bulges are resolved and have typical ranges of 0.5-2 k pc while the 



mairi dust heating by AGN (500-2000K) is generated on scales of (10-100) pc fjSoifer et al. 



20011 ). Furthermore, the central light excess expected from an AGN contributes on average 



less than 5% to the total bulge luminosity for our sample (see Fig. [18]) . 
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5. Summary 

We have studied the nuclear stellar properties of 73 LIRG systems (\og[L iji/LQ] = 
11.4 — 12.5) using high-resolution H-band images obtained with HST NICMOS revealing 
a large variety of double and triple nuclei merger systems. To investigate the effect of the 
merging process on the evolution of galaxies, such as the growth of the central BH or stellar 
concentration, we studied several morphological parameters as a function of merger stage 
using GALFIT. The main results are summarized as follows: 

1. The fraction of LIRG systems with at least two interacting nuclei is 63%. Triple nu- 
clei systems are found for at least 3 (possibly 5) LIRG systems. The comparison of 
our NICMOS images with comparable resolution HST B-band images revealed that 
roughly 50% of the double nuclei are obscured by dust, with 10 of the obscured nu- 
clei responsible for the primary mid-IR emission (MIPS 24/im). This implies strong 
limitations on the ability to detect the true nuclear structures of luminous infrared 



the LIRG population at local and high-re dshift (e.g 



Kim et al. 


2002; 


Bell et al. 


2005 



Melbourne et all 120051 : lOasvra et al.ll2008f ). 



2. ULIRGs (log[L//j;/LQ] > 12.0) have significantly smaller nuclear separations than LIRGs 
(log[L7;^/LQ]= 11.4 — 12.0) with a median value of 1.2 kpc and 6.7 kpc, respectively. 
In our sample, merger (regardless of whether LIRG or ULIRG) seem to be prevalent 
at two time scales (based on the projected nuclear separation and mass ratio of the 
nuclei): First, at a remaining merger time scale of 0.3 < [t — tmerg] < 1-3 Gyr (53% of 
mergers in our sample), and second, at [t — tmerg] ~ (26%), likely representing the 
first passage of interacting galaxies and the final nuclear coalescence, respectively. An 
average post-merger time (LIRG phase after the nuclei merged) of 0.3 Gyr is estimated. 

3. The bulge luminosity shows, on average, a small increase (about a factor of 1.8) towards 
late merger stages. Interestingly, the LIRGs in our sample that show no interaction 
features have on average a significantly larger bulge luminosity (a factor of two) than 
interacting LIRGs. This is likely to be a selection effect for more intrinsically luminous 
non-merging galaxies among the LIRG sample. 

4. The bulge luminosity surface density L Bulge / R%uige increases significantly along the 
merger sequence which is primarily due to a decrease of the bulge radius towards 
late merger stages. These findings are in agreement with models that include gas 
inflow which can dramatically affect the measured bulge properties during the merging 
process. No significant dependence of the bulge Sersic index is found as a function of 
merger stage. 
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5. Although no significant correlation between the 6.2//m PAH EQW and the core excess 
light in the H band is found, LIRGs with the largest core excess exhibit slightly smaller 
PAH EQWs, suggesting more hot dust emission and hence AGN activity. 

6. At least seven LIRGs in our sample exhibit warm dust emission (> 70% of its total 
infrared luminosity) that is clearly offset form the nucleus. Instead, their emission 
regions seem to be associated with spiral arms, a possible secondary nucleus, or the 
region between the merging nuclei. 

We would like to thank Philip F. Hopkins for valuable discussions of our results in 
terms of model predictions and numerical simulations. This research has made use of the 
NASA/IPAC Extragalactic Database (NED) and Infrared Science Archive which is operated 
by the Jet Propulsion Laboratory, California Institute of Technology, under contract with 
the National Aeronautics and Space Administration. Support for this work was provided 
through grant HST GO 11235.01-A by NASA from the Space Telescope Science Institute, 
which is operated by the Association of Universities for Research in Astronomy, Inc., under 
NASA contract NAS 5-26555. 
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Table 1. Sample Overview 
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85 
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01 


44 
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64 
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51 
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38 
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28 


87.9 


9726 


NGC6090 


16 


11 


40 


60 


+ 52:27:25.0 


11 


58 


137 


7219 


2MASXJ16191179-0754026 


16 


19 


11 


80 


-07:54:03.0 


11 


62 


128 


11235 


ESO069-IG006N 


16 


38 


11 


84 


-68:26:10.4 


11 


98 


212 


11236 


ESO069-IG006S 


16 


38 


13 


48 


-68:27:19.3 


11 


98 


212 


11236 


IRAS 16399-0937 


16 


42 


40 


20 


-09:43:14.0 


11 


63 


114 


11236 


NGC6240 


16 


62 


68 


70 


+02:24:04.0 


11 


93 


116 


7219 


IRASF17132+5313 


17 


14 


20 


24 


+53:10:30.8 


11 


96 


232 


11235 


IRAS-F17138-1()17 


17 


16 


35 


60 


-10:20:38.0 


11 


49 


84 


10169 


IRAS-F17207-0014 


17 


23 


22 


20 


-00:17:02.0 


12 


46 


198 


7219 


IRAS 18090+0130 


18 


11 


38 


42 


+01:31:38.8 


11 


65 


134 


11236 


IC4689 


18 


13 


40 


28 


-57:44:53.5 


11 


62 


81.9 


11236 


IRAS 18293-3413 


18 


32 


41 


22 


-34:11:26.7 


11 


88 


86 


11236 


NGC6670W 


18 


33 


33 


77 


+69:63:16.6 


11 


65 


129.5 


11236 
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Table 1 — Continued 



Name 


RA 


DEC 




D 


HST ID 




(J2000) 


(J2000) 


[Mpc] 




NGC6670E 


18:33:37.70 


+59:53:23.0 


11.65 


129.5 


11235 


NGC6786S 


19:10:53.90 


+ 73:24:37.0 


11.49 


113 


11235 


UGC11415 


19:10:53.90 


+73:24:37.7 


11.49 


113 


11235 


NGC6786N 


19:11:04.44 


+73:26:37.4 


11.49 


113 


11236 


ESO693-IG008 


19:14:30.98 


-21:19:06.8 


11.93 


222 


11236 


IRAS-F19297-0406 


19:32:22.28 


-04:00:01.1 


12.45 


395 


7896 


IRAS19542+1110 


19:66:35.39 


+ 11:19:03.0 


12.12 


295 


11235 


IRAS20351+2521 


20:37:17.80 


+25:31:38.0 


11.61 


151 


11235 


IIZW()9B 


20:57:23.94 


+ 17:07:39.5 


11.94 


161 


11235 


BSO286-IG019 


20:58:26.80 


-42:39:00.0 


12.06 


193 


11235 


IRAS21101+5810 


21:11:29.82 


+68:23:07.2 


11.81 


174 


11236 


BSO239-IG002 


22:49:39.90 


-48:50:58.0 


11.84 


191 


11236 


IRAS-F22491-1808 


22:51:49.30 


-17:52:24.0 


12.2 


351 


7219 


NGC7469 / IC5283 


23:03:17.88 


+08:53:39.3 


11.65 


70.8 


11235 


ESO148-IG002 


23:15:46.79 


-59:03:13.0 


12.06 


199 


7896 


IC5298 


23:16:00.71 


+ 25:33:24.0 


11.6 


119 


11235 


BSO077-IG014 


23:21:04.44 


-69:12:54.8 


11.76 


186 


11235 


NGC7674 


23:27:66.70 


+08:46:46.0 


11.56 


125 


11236 


IRASF23366+3604 


23:39:01.30 


+36:21:09.8 


12.2 


287 


11235 


IRAS23436+5257 


23:46:05.59 


+ 53:14:01.0 


11.57 


149 


11235 


UGC12812W / MRK0331 


23:51:18.73 


+ 20:34:42.9 


11.5 


79.3 


11235 


UGC12812E / MRK0331 


23:51:26.80 


+20:35:10.0 


11.5 


79.3 


11235 



Note. — Summary of the properties of our 73 LIRG systems (88 pointings). Column (1): Source Name from NED, Column (2): right ascension 
(J2000), Column (3): source declination (J2000), Column (4): The luminosity distance in Mpc (adopting Hq = 70km s~^ Mpc), as provided by 
NED, Column (5): The total infrared luminosity in logio Solar units, Column (6): The data origin given by the ID of the observational program. 



Table 2. Bulge Parameters and Merger Classification 



Name 




DECj,,,,,,, 


M, Class. 


Lb 










LpsF 


Ci.VLFIT romp. 




(.1201)1)) 


(.T2(H)(l) 




l.)gi().L_.,J 


[klK- 






ll()^"xL_.,,.ki.c--J 







NGC0034 


00: 


ill: 


:06. 


.544 


-12:06:27.24 


6 


11, 


.12 ± 10.08 


0, 


.13 ± 0.0 


1. 


,87 


± 0.01 


8.66 ± 0.09 


777.34 


8.97 


2 


ARP256S 


00; 


:18: 


:50, 


.870 


-10:22:36.46 


3 


10, 


.87 ± 9.54 


0, 


.35 ± 0.03 


2, 


,62 


± 0.16 


8.27 ± 0.07 


59.76 


9.35 


1 


MCG-03-04-014 


01 


10 


08 


912 


-16:51:09.58 





11 


,32 ± 9.98 


1 


,05 ± 0.04 


2 


07 


± 0.04 


8.79 ± 0.1 


18.85 


9.44 


2 


CGCG436-030 


01 


20 


()2 


632 


+ 14:21:42.26 


2 


10 


,84 ± 9.4 


0, 


,22 ± 0.02 


6 


76 


± 0.42 


8.23 ± 0.07 


136.95 


9.46 


2 


2MASX.T01385289-1027113 


01 


:38: 


:52, 


.857 


-10:27:11.38 


5 


10, 


,42 ± 8.39 


0, 


,38 ± 0.0 


1 


39 


± 0.02 


7.75 ± 0.08 


18.72 


8.66 


2 


IIIZW03,S 


01: 


:44: 


:30, 


,546 


+ 17:06:09.04 


3 


11 


.15 ± 9.11 


2 


.07 ± 0.02 


2. 


32 


± 0.04 


8.59 ± 0.08 


3.29 


8.66 


2 


NGC0695 


01: 


:51: 


:14, 


.333 


+22:34:56.02 





11 


.41 ± 10.07 


1 


.85 ± 0.11 


3, 


,73 


± 0.07 


8.9 ± 0.1 


7.54 


8.91 


2 


MRK1034W 


02: 


:23: 


:18. 


.963 


+32:11:18.69 


2 


11 


.1 ± 9.54 


2, 


.08 ± 0.1 


7. 


.75 


± 0.11 


8.64 ± 0.08 


2.93 


8.82 


3 


MRK1034E 


02: 


:23: 


:21 


.949 


+32:11:48.83 


2 


11 


.56 ± 9.52 


2, 


.91 ± 0.02 


1, 


,14 


± 0.01 


9.07 ± 0.12 


4.29 


_ 


3 


IR.ASF033.59+1523 


03: 


:38: 


:46, 


,351 


+ 15:32:54.69 


3 


10 


,99 ± 8.96 


2, 


,82 ± 0.07 


7, 


45 


± 0.07 


8.41 ± 0.07 


1.24 


_ 


2 


ESO,S,S0-IG02N 


04: 


:21 


:19, 


,977 


-18:48:39.39 


2 


10 


,59 ± 8.55 


0, 


,19 ± 0.0 


1 


37 


± 0.01 


7.94 ± 0.07 


107.86 


8.88 


2 


ESO.S,S0-IG02S 


04: 


:21 


:20, 


,017 


-18:48:57.17 


2 


10, 


,5 ± 9.46 


0, 


,2 ± 0.0 


1 


29 


± 0.01 


7.84 ± 0.09 


81.35 


_ 


2 


NGC1614 


04: 


:34: 


:00, 


,015 


-08:34:45.13 


5 


10, 


,9 ± 8.86 


0, 


,08 ± 0.0 


1. 


55 


± 0.02 


8.31 ± 0.07 


1309.26 


10.23 


3 


BSO203-IG001 


04: 


:46: 


:48, 


.986 


-48:33:35.10 


3 


10, 


.45 ± 8.42 


2, 


.41 ± 0.05 


2, 


,94 


± 0.04 


7.79 ± 0.08 


0.49 


_ 


1 


ESO203-IG001 


04: 


:46: 


:49, 


.537 


-48:33:29.90 


3 


10, 


.87 ± 8.84 


1 


.37 ± 0.02 


2. 


,02 


± 0.01 


8.27 ± 0.07 


3.93 


8.3 


2 


VII-Zw-031 


05: 


:16: 


:46, 


.502 


+79:40:12.85 





11 


.34 ± 10.38 


1 


.09 ± 0.18 


4, 


,47 


± 0.28 


8.81 ± 0.11 


18.35 


_ 


2 


ESO2,5,5-IG007N 


06: 


:27: 


:21 


,754 


-47:10:35.78 


3 


11 


,25 ± 10.36 


0, 


,61 ± 0.1 


1, 


,35 


± 0.2 


8.71 ± 0.11 


47.82 


_ 


1 


ESO2,S,S-IG007N 


06: 


:27: 


:22 


,587 


-47:10:46.79 


3 


10 


,84 ± 9.76 


0, 


,96 ± 0.09 


1 


6 ± 0.13 


8.24 ± 0.08 


7.57 




1 


ESO2,S,'3-IG007S 


06: 


:27: 


:23, 


,090 


-47:11:02.73 


3 


11 


,11 ± 9.38 


9, 


,37 ± 0.6 


3. 


12 


± 0.1 


8.55 ± 0.08 


0.15 




2 


AM0702-601N 


07: 


:03: 


:24, 


,257 


-60:15:22.53 


1 


10, 


,71 ± 9.93 


1 


,81 ± 0.38 


4, 


,45 


± 0.39 


8.09 ± 0.11 


1.58 




3 


AM0702-601S 


07: 


:03: 


:28, 


.627 


-60:16:44.67 


1 


11 


.43 ± 10.4 


2, 


.14 ± 0.03 


3, 


,14 


± 0.04 


8.92 ± 0.11 


5.92 




2 


NGC2623 


08: 


:38: 


:24, 


.094 


+25:45:16.75 


6 


10, 


.61 ± 8.57 


0, 


.13 ± 0.0 


3. 


,39 


± 0.04 


7.97 ± 0.07 


225.78 


8.45 


2 


ESO060-IG016 


08: 


:52: 


:32, 


.131 


-69:01:55.50 


3 


10, 


.84 ± 9.88 


1 


.27 ± 0.19 


2. 


,04 


± 0.22 


8.24 ± 0.09 


4.33 




1 


IRAS-F08572+3915 


09: 


:00: 


:25, 


,593 


+39:03:50.60 


3 


9.92 ± 7.88 


0, 


,44 ± 0.0 


0, 


,91 


± 0.02 


7.16 ± 0.13 


4.18 




2 


UGC05101 


09: 


:35: 


:51, 


,629 


+ 61:21:11.89 


5 


11 


,48 ± 9.75 


0, 


,84 ± 0.02 


3, 


,05 


± 0.07 


8.98 ± 0.11 


43.14 


10.12 


2 


IRAS-F10565+2448 


10: 


:59: 


:18, 


,146 


+ 24:32:34.54 


2 


11 


,65 ± 10.61 


1 


,0 ± 0.01 


2 


73 


± 0.02 


9.18 ± 0.13 


44.7 


10.09 


1 


IRAS-F11231+1466 


11 


:25: 


:45, 


,077 


+ 14:40:36.07 


1 


11 


,4 ± 9.84 


3, 


,29 ± 0.15 


.5 


06 


± 0.08 


8.88 ± 0.1 


2.31 




2 


NGC3690E 


11 


:28: 


:33, 


.690 


+58:33:46.35 


3 


10, 


.85 ± 9.82 


0, 


.68 ± 0.07 


3, 


,54 


± 0.13 


8.25 ± 0.08 


15.57 


_ 


2 


IRAS-F12112+0306 


12: 


:13: 


:46, 


.904 


+02:48:39.26 


4 


11 


,0 ± 9.44 


0, 


.32 ± 0.02 


5. 


.98 


± 0.4 


8.42 ± 0.08 


99.63 


8.91 


1 


WKK0787 


12: 


:14: 


:22, 


.098 


-56:32:33.26 





11 


.43 ± 10.4 


1 


.44 ± 0.01 


2, 


,95 


± 0.02 


8.93 ± 0.11 


13.07 


9.6 


1 


WKK2031 


13: 


:15: 


:06, 


.339 


-55:09:22.72 


5 


11 


,56 ± 9.22 


0, 


.46 ± 0.0 


1, 


11 


± 0.01 


9.07 ± 0.11 


170.09 


10.43 


2 


UGC08335W 


13: 


:15: 


:30, 


.799 


+62:07:45.31 


2 


10, 


,43 ± 9.3 


0, 


,49 ± 0.06 


4, 


72 


± 0.34 


7.76 ± 0.09 


11.31 




2 


UGC08335E 


13: 


:15 


:3,5, 


,023 


+ 62:07:28.86 


2 


11 


,23 ± 9.19 


0, 


,76 ± 0.01 


3 


72 


± 0.04 


8.69 ± 0.09 


29.4 


8.44 


2 


UGC08387 


13: 


:20 


:3,''j 


,319 


+ 34:08:22.39 


4 


10 


,51 ± 8.47 





,55 ± 0,01 


1 


27 


± 0.02 


7.85 ± 0.08 


10.55 




2 


NGC5256 


13: 


:38: 


:17, 


.270 


+ 48:16:32.10 


3 


10, 


.87 ± 8.84 


0, 


.62 ± 0.0 


1. 


,02 


± 0.01 


8.28 ± 0.07 


19.45 




2 


NGC5256 


13: 


:38: 


:17, 


.764 


+48:16:41.19 


3 


10, 


.96 ± 9.77 


0, 


.76 ± 0.07 


3. 


,47 


± 0.12 


8.37 ± 0.08 


15.65 


8.39 


2 


NGC5257 


13: 


:39: 


:53, 


.767 


+00:50:26.78 


2 


11 


.03 ± 9.6 


2, 


.29 ± 0.16 


6, 


,22 


± 0.16 


8.46 ± 0.08 


2.06 




2 


IRAS-F14348-1447 


14: 


:37: 


:38, 


,285 


-15:00:24.05 


4 


10, 


.92 ± 9.48 


0, 


,36 ± 0.01 


1, 


51 


± 0.06 


8.32 ± 0.07 


63.44 


9.64 


2 


IRAS-F14348-1447 


14: 


:37: 


:38, 


,402 


-15:00:21.13 


4 


11 


,57 ± 10.13 


5 


,87 ± 0.36 


4 


7 ± 0.12 


9.08 ± 0.12 


1.08 


8.74 


1 


IRAS-F14378-3651 


14: 


:40: 


:59, 


.013 


-37:04:31.89 


6 


11 


.19 ± 10.16 


0, 


.24 ± 0.0 


1, 


,43 


± 0.01 


8.65 ± 0.1 


262.25 




2 


UGC09618S 


14: 


:67: 


:00. 


.327 


+24:36:24.02 


1 


11 


.66 ± 9.53 


7, 


.23 ± 0.06 


3. 


,39 


± 0.01 


9.08 ± 0.12 


0.7 


8.95 


2 


ESO099-G004 


16: 


:24: 


:67, 


.964 


-63:07:29.67 


3 


10, 


.92 ± 8.89 


0, 


.4 ± 0.01 


5. 


,42 


± 0.16 


8.33 ± 0.07 


51.13 




1 


IRAS-F15250+3608 


16: 


:26: 


:69, 


.423 


+35:58:37.22 


6 


10, 


.86 ± 9.86 


0, 


.66 ± 0.08 


3. 


,21 


± 0.24 


8.26 ± 0.09 


22.81 




2 


NGC6090 


16: 


ill 


:40, 


.918 


+52:27:27.22 


4 


10, 


.07 ± 8.64 


0, 


.16 ± 0.01 


2. 


,16 


± 0.12 


7.34 ± 0.11 


46.12 




2 


2MASX.I16191179-07,54026 


16: 


:19: 


:11 


,787 


-07:54:03.02 


5 


11 


.33 ± 9.6 


4, 


.98 ± 0.17 


4, 


,09 


± 0.08 


8.81 ± 0.1 


0.86 




2 


ESO069-IG00BS 


16: 


:38: 


:13, 


,472 


-68:27:16.83 


2 


11 


,25 ± 9.52 


1 


,31 ± 0.03 


3, 


24 


± 0.03 


8.72 ± 0.09 


10.48 




2 


IRAS16399-0937 


16: 


:42: 


:40, 


.141 


-09:43:13.20 


3 


10, 


.71 ± 9.16 


1 


.03 ± 0.02 


1. 


,26 


± 0.07 


8.09 ± 0.07 


4.91 


8.69 


2 


IRAS16399-0937 


16: 


:42: 


:40, 


.178 


-09:43:18.74 


3 


10, 


.63 ± 9.55 


0, 


.61 ± 0.06 


2. 


,93 


± 0.15 


7.99 ± 0.08 


16.38 




2 


NGC6240 


16: 


:62: 


:68, 


.886 


+02:24:03.20 


4 


11 


.29 ± 10.25 


0, 


.23 ± 0.0 


1. 


,6 ± 0.01 


8.76 ± 0.1 


376.66 


9.45 


1 


NGC6240 


16: 


:62: 


:68, 


.934 


+02:24:04.92 


4 


10, 


.81 ± 9.38 


0, 


.2 ± 0.01 


1. 


,82 ± 0.07 


8.2 ± 0.07 


155.08 




2 



Table 2 — Continued 



N anic 


RA£ 








iVi . v.^iass. 






Rb 












^PSF 


'.jrt.ijfil coiiip. 






(J2000) 


J ^UUU J 




1 n tt 1 

10glU[LQj 


[kpc] 










10glO[AiQj 


[lU Xij0/Kpc J 


1 „„i n n 1 

10glU[ijQj 




IRASFl 7132+5313 


17 


■ 14' 


19 


802 


+53:10:28.89 


2 


11.09 ± 9.36 


Q 


71 it 02 


5 


.41 


± 


0.08 


8.53 ± 0.08 


24.46 


9.01 


2 


IRASFl 7132+5313 


17 


14 


20 


453 


+53:10:31.99 


2 


11.04 ± 10.4 




16 it 31 


2 


.81 


± 


0.32 


8.47 ± 0.14 


8.12 




3 


IRAS-F 17207-0014 


17 


■23' 


21 


953 


-00:17:00.75 


5 


10.39 ± 9.06 


Q 


.oo HI u.ux 





.96 


± 


0.05 


7.72 ± 0.09 


23.41 




2 


IRAS 18090+0130 


18 


: 11; 


38. 


A12 


+01:31:39.99 


2 


10.54 ± 8.5 


Q 


.26 it 0.0 


1 


.12 


± 


0.03 


7.89 ± 0.08 


51.38 




2 


IC4689 


18 


:13: 


40. 


386 


-57:44:54.18 


2 


10.32 ± 8.29 


Q 


.22 it 0.0 


2 


.86 


± 


0.03 


7.64 ± 0.09 


45.48 


8.15 


2 


IRAS18293-3413 


18 


■32' 


41 


139 


-34:11:27.61 


1 


11.58 ± 9.84 




TQ A- n OS 


5 


.67 


± 


0.13 


9.09 ± 0.12 


11.86 




2 


NGC6670W 


18 


:33; 


;33. 


.972 


+59:53:17.49 


2 


11.09 ± 9.05 


2 


.04 ± 0.02 


1 


.78 


± 


0.02 


8.52 ± 0.08 


2.94 




2 


NGC6670E 


18 


:33; 


;37. 


721 


+59:53:23.07 


2 


10.65 ± 9.62 





.88 ± 0.0 


0. 


.82 


± 


0.01 


8.02 ± 0.09 


5.88 




2 


NGC6786N 


19 


:11; 


;04. 


.302 


+73:25:33.42 


2 


10.59 ± 9.56 





.1 ± 0.0 


2 


.87 


± 


0.05 


7.95 ± 0.09 


398.67 


9.14 


2 


ESO593-IG008 


19 


:14: 


;31 


.099 


-21:19:09.33 


4 


11.35 ± 10.09 


1 


.28 ± 0.1 


2 


.37 


± 


0.14 


8.83 ± 0.1 


13.61 


7.49 


1 


IRAS19542+1110 


19 


:56: 


;35. 


787 


+ 11:19:05.09 





11.5 ± 10.17 


2 


.04 ± 0.06 


2 


.95 


± 


0.11 


9.01 ± 0.11 


7.7 


11.56 


2 


IIZW096 


20 


:57: 


24. 


475 


+ 17:07:39.89 


3 


10.09 ± 9.38 





.59 ± 0.16 


3. 


.58 


± 


0.44 


7.37 ± 0.15 


3.59 




2 


ESO286-IG019 


20 


:58: 


;26. 


.801 


-42:39:00.20 


5 


11.18 ± 9.15 





.7 ± 0.01 


1 


.92 


± 


0.02 


8.64 ± 0.08 


31.2 


9.04 


2 


IRAS21101+5810 


21 


:11: 


;30. 


.392 


+58:23:03.38 


2 


10.42 ± 8.69 





.39 ± 0.01 


2 


.16 


± 


0.02 


7.75 ± 0.08 


17.5 




2 


ESO239-IG002 


22 


:49: 


;39. 


.889 


-48:50:58.20 


5 


11.18 ± 10.15 





.16 ± 0.0 


3. 


.02 


± 


0.02 


8.63 ± 0.1 


571.59 


9.96 


2 


IRAS-F22491-1808 


22 


:51: 


;49. 


.354 


-17:52:24.08 


4 


10.14 ± 8.8 





.32 ± 0.02 


0. 


.97 


± 


0.15 


7.42 ± 0.11 


13.39 




1 


NGC7469 / IC5283 


23 


:03: 


17. 


.964 


+08:53:36.72 


2 


9.78 ± 8.04 





.24 ± 0.01 


1 


.92 


± 


0.05 


7.0 ± 0.14 


10.42 




3 


IC5298 


23 


:16; 


;00. 


.682 


+25:33:24.02 





11.16 ± 9.42 


1 


.17 ± 0.01 


1 


.17 


± 


0.03 


8.6 ± 0.08 


10.41 


9.63 


3 


ESO077-IG014 


23 


:21: 


;03. 


.620 


-69:13:00.99 


2 


11.32 ± 9.28 


3 


.44 ± 0.03 


2 


.33 


± 


0.02 


8.79 ± 0.09 


1.75 




3 


ESO077-IG014 


23 


:21: 


;05. 


.329 


-69:12:47.30 


2 


11.47 ± 9.91 


2 


.55 ± 0.11 


4. 


.68 


± 


0.08 


8.96 ± 0.11 


4.52 


9.25 


3 


NGC7674 


23 


:27: 


;56. 


717 


+08:46:44.40 


2 


11.07 ± 10.03 





.37 ± 0.0 


4. 


.27 


± 


0.04 


8.5 ± 0.09 


86.46 


10.37 


3 


IRASF23365+3604 


23 


:39: 


;01 


.269 


+36:21:08.54 


5 


10.96 ± 9.92 





.26 ± 0.0 


1 


.65 


± 


0.02 


8.37 ± 0.09 


134.53 




2 


IRAS23436+5257 


23 


:46: 


;05. 


512 


+53:14:01.40 


4 


10.58 ± 8.85 





.24 ± 0.0 


1 


.5 ± 0.03 


7.94 ± 0.08 


67.86 


9.66 


2 


UGC12812W 


23 


:51: 


18. 


.675 


+20:34:41.59 


1 


9.57 ± 8.86 


1 


.64 ± 0.17 


1 


.5 ± 0.03 


6.77 ± 0.19 


0.14 




2 


UGC12812E 


23 


:51: 


;26. 


741 


+20:35:10.41 


1 


11.06 ± 9.5 


2 


.3 ± 0.15 


7. 


.43 ± 0.19 


8.5 ± 0.08 


2.19 




3 



Note. — Column (1): Source Name from NED, Column (2): right ascension (J2000) of bulge center, Column (3): declination {J2000) of bulge center, Column (4): The merger 
classification as described in ^ |3.2| Column (5): The apparent bulge luminosity with fit error in mag , Column (6): The bulge radius with fit error in kpc, Column (7): The Sersic index of 
the bulge radial profile with fit error, Column (8): The mass of the central black hole estimated from the bulge luminosity (based on the relation of lMarconi Huntir2003l ). Column (9): 
The bulge surface density ^Bulcje/^^i^igQ2 i Column (10): The luminosity of the central PSF, Column (11): The number of GALFIT components (excluding PSF and background). 



Table 3. Merger Stage Results 



Bulge property 


Statistics 


1 


2 


3 


4 


5/6 


IL log[L0l 


Mean 


11.38 


11.09 


10.88 


11.06 


11.14 


Median 


11.41 


11.04 


10.87 


10.87 


11.12 




Mean Error 


9.95 


9.75 


9.59 


9.73 


9.74 




Std. Dev. 


9.87 


9.97 


9.72 


9.8 


9.71 


R [kpc] 


Mean 


3.09 


1.08 


1.54 


0.95 


0.71 




Median 


2.22 


0.74 


0.76 


0.32 


0.33 




Mean Error 


0.14 


0.04 


0.09 


0.05 


0.02 




Std. Dev. 


1.92 


0.99 


2.08 


1.67 


1.25 


Ha 


Mean 


4.86 


3.26 


2.8 


2.39 


2.2 




Median 


4.76 


2.84 


2.62 


1.71 


1.87 




Mean Error 


0.14 


0.09 


0.12 


0.11 


0.05 




Std. Dev. 


1.45 


1.93 


1.61 


1.55 


0.97 


LSD log[lOl°XL0/kpc^] 


Mean 


10.61 


11.69 


11.18 


11.93 


12.44 


Median 


10.35 


11.16 


10.69 


11.74 


12.13 




Mean Error 


9.69 


10.64 


10.43 


10.89 


11.2 




Std. Dev. 


10.58 


11.93 


11.27 


12.03 


12.57 



Note. — The mean, median, mean systematical error (from GALFIT), and the standard deviation of the bulge properties (NIR luminosity 
L, radius R, Scrsic index n^, and luminosity surface density LSD) for all merger stages. Note that we have combined merger stage 5 and 6 as 
post-merger stage LIRGs. 



Table 4. KS and MWU Test Results 







■ly Tc^l 


1 - 2 


1 - :i 


i - 4i 


1 -r,:ii 2 




/(i :! - 4 :l - 4 - r,/ii 



L 


K-S p 




0.01 


0.01 


0.02 


0.04 


0.11 


0.74 


0.54 


0.54 


0.11 


0.74 




MWU 


P 


0.01 


0.01 


0.01 


0.03 


0.48 


0.68 


0.43 


0.8 


0.09 


0.42 


R 


K-S p 




0.01 


0.01 


0.01 


0.01 


0.15 


0.21 


0.12 


0.01 


0.01 


0.88 




MWU 


P 


0.01 


0.01 


0.01 


0.01 


0.31 


0.25 


0.09 


0.01 


0.01 


0.85 


na 


K-S p 




0.09 


0.05 


0.01 


0.01 


0.54 


0.19 


0.32 


0.42 


0.36 


0.96 




MWU 


P 


0.19 


0.05 


0.04 


0.01 


0.55 


0.19 


0.17 


0.38 


0.41 


0.95 


LSD 


K-S p 




0.06 


0.73 


0.01 


0.01 


0.39 


0.28 


0.04 


0.06 


0.01 


0.37 




MWU 


P 


0.06 


0.76 


0.01 


0.01 


0.09 


0.18 


0.02 


0.02 


0.01 


0.21 



Note. — The Kolmogorov-Smirnov (KS) and Mann- Whitney U (MWU) test results of the bulge properties (NIR luminosity L, radius R, Scrsic 
index as, and luminosity surface density LSD) for all combinations of merger stage populations. Note that we have combined merger stage 5 and 
6 as post-merger stage LIRGs. The p values listed are the probability levels of the hypothesis that the two samples have the same distribution. 
Large p values indicate that the two samples come from the same sample distribution while small p-values (p;^ 0.15, underlined) are not drawn 
from the same parent population. The KS test tests differences in the shapes of the distributions of the two groups, not just the locations of the 
ranks as in the Mann- Whitney U test. 
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